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Abstract 

The double lepton pair decay modes of the Kl meson are analyzed including 
all contributions of order p 6 in Chiral Perturbation Theory. The experimen- 
tally established e + e~e + e~ mode and the recently observed e + e~ fi + mode 
are discussed in detail. 
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In this note we reconsider the double lepton pair decays Kl — > £ + £~£' + £'~ that were 
studied many years ago by Miyazaki and Takasugi The two modes of main interest are 
the Kl — > e + e~e + e~ and the Kl — > e + e~fi + fi~ decays. The first proceeds at a rate well 
determined experimentally [Q] with a branching ratio equal to 4.1±0.8 x 10~ 8 , and its interest 
resides partly in the possibility of determining the form factor associated with the virtual 7 
couplings to the Kl as well as the interesting interference effect derived from the identical 
leptons in the final state. Both issues are analyzed here, and the main conclusions are that 
the form factor effect is small (4%), requiring a substantial experimental improvement over 
the current error of about 20%, and the interference effect is tiny (0.5%) and most likely 
beyond experimental access. The second decay has been recently observed for the first time 
by the E799 collaboration at Fermilab ||. This experiment quotes a branching ratio of a 
few parts per billion, a remarkable improvement over the previous upper bound of 4.9 x 10 -6 
0. We discuss in this case the effects of a non-trivial form factor, that according to our 
result produces an increase of the branching ratio by 30%, which could therefore be tested 
in the future as the experimental accuracy in this mode improves. 

In the following we neglect a contribution to the amplitude due to CP violation. This 
contribution is of order p 4 but suppressed by the CP violating mixing parameter e which is 
of order 10~ 3 (Re(e) = 1.6 x 10~ 3 , e ~ 43.5°). Thus, the CP violating contributions to the 
widths considered here are suppressed by a factor approximately equal to | e | A 2 x jM\ ~ 1%, 
where A x = 4tt F n . Thus, in the limit of CP conservation, the Kl — > 7*7* amplitude has 
the most general form: 

A(K L - 7*7'*) = Ht, ^e' P K- (1) 

Here 7* and 7'* are virtual photons with respective invariant mass squared t = k 2 and 
t' = k' 2 . !F(t,t') is the form factor of order p 6 in the chiral expansion that has a t and t' 
independent contribution Fx, due to the 7r°, 77 and rj' poles ||, plus a t and t' dependent 
contribution i^t, t') from one chiral loop plus counterterms. F\ is entirely fixed up to its 
sign by the Kl — > 77 decay, and F 2 (t 1 t') is given by ||: 
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F 2 (t, if) = {-(a a + 2a 4 ) D(t, f , //) + C{fi (t + t')} , (2) 

where the counterterm has eliminated an UV divergence proportional to (t + t'), F n = 93 
MeV is the pion decay constant, and 

£>(t, f, /i) = (t + [- - (log -f + log -*)] 
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F(M^ i) + F[M K , t) + F(M*, + F(Mk, t') , (3) 



with the chiral logarithms contained in the function F(m 2 , t): 



F(m , i) = (1 - -) J log-^— , - 2 ra , 

( 4 ) 

The coefficient C(/j,) in the counterterm is determined by the fit to the Dalitz decays as 
discussed in ||. There are two scenarios distinguished by the sign of Fi, where the relative 
size of the counterterm to the chiral logarithms is different. In one scenario the counterterm 
(defined at the subtraction scale /i = M p ) provides more than 90% of the contribution to 
the form factor's slope, while in the other that fraction is reduced to about 75%. 

Since ^-"(0, t) = J-(t, 0) = F(t), where F(t) is the form factor of the Dalitz decays studied 
in 0, we use the results obtained in that reference by fitting the data M: 



Scenario 1 : 



F l = 0.89 — -, a 2 + 2a 4 = -0.3 ± 0.3, C(p = M p ) = 14.2 ± 7.3 
Scenario 2 : 

c 

F x = -0.89 " em ^ 8 , a 2 + 2a 4 = 1.5 ± 0.3, C(jjt = M p ) = -10.3 ± 7.3 (5) 

ZTlFn 

Here, C% = 3.12 x 10 -7 is the octet coupling in the non-leptonic weak interaction effective 
Lagrangian of order p 2 . 

The decay amplitude has one piece if the final lepton pairs are different, and two pieces 
if they are identical. In this latter case the two amplitudes are: 
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Ai = e 2 F(t, t') <T<" (g±±^±±^ 

x u(P-) 7^ u(p+) u(p'_) 7 P i>(p' + ) . (6) 

ss' 

x 7^ t>(p + ) tZ(p_) 7p v(p+) , (7) 

where p + is the momentum of £ + , etc., and 

t = ( P+ +p~) 2 , t' = {p' + +p'_)\ 

s = (p + +p'_) 2 , s '=(p' + +p^ 2 . (8) 

In the case of distinguishable lepton pairs only one appears, say A\. The decay width is 
obtained by summing over the lepton spins and integrating over the four particle phase 
space. We checked the results obtained in |J and refer to it for further details. 
In the Kl — * £ + £~£ + £~ decay we have: 

T(K L -> £+£-£+£-) = ~ T 1 (K L -> £ + £~£ + £~) + i T 2 (K L -> £+£-£+£-) (9) 

with 

Tl = / E Ml or 2 T , 

spins 

r 2 = / Re(E A^)^, (10) 

spins 

where c?$ represents the four body phase space volume element. On the other hand, in the 
Kl — > e + e~fi + fi~ decay we have instead: 

T(K L -> e+e-/iV) = -> e+e>V) (H) 

For convenience we refer to the — * 77 width. Defining p = T/Y{Ki — > 77) and 
^interference = | Tz/V^Kl — > 77), we obtain the results shown in the table: 
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Decay mode 


Analysis 


Pinterference 


P 




Ref i 


-0.35 x KT 5 


5.89 x 10~ 5 




No form factor 


-0.036 x 10~ 5 


6.26 x lO" 5 




With form factor 


-0.048 x 10~ 5 


6.50 x 10~ 5 






-0.047 x 10~ 5 


6.48 x 10~ 5 




Ref i 





1.42 x 10~ 6 




No form factor 





1.71 x 10~ 6 




With form factor 





(2.20 ±0.25) x 10~ 6 









(2.18 ±0.25) x 10~ 6 




Ref i 


-0.051 x 10~ 9 


0.946 x 10" 9 




No form factor 


-0.051 x 10~ 9 


0.93 x 10~ 9 




With Form Factor 


-0.077 x 10~ 9 


(1.30 ±0.15) x 10~ 9 






-0.072 x 10~ 9 


(1.35 ±0.15) x 10- 9 



TABLE: The results of Ref [T|] correspond to a point like form 
factor. The results of this work with a form factor are given 
respectively for the two scenarios of equation (5). 



In the e + e~e + e~ mode there is a small effect due to the form factor leading to an increase 
of the width by about 4%. Given the current experimental error of almost 20%, it seems 
that a test of the form factor can be achieved in the foreseeable future. On the other hand, 
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we find that the interference term due to the identity of particles is a factor ten smaller 
than that obtained before and it represents a correction of 0.5%, which seems beyond 
experimental access. Our prediction is consistent with the experimental rate: BR(Kl — > 
e + e-e + e-)| T heor y = 3.85 x 1(T 8 , and BR(K L -> e+ e" e+ e-)| Exp . = (4.1 ± 0.8) x 10~ 8 . The 
relative size of the interference effect is larger in the fi + fi~{i + {i~ mode, but alas, the total 
branching ratio for this decay is predicted to be about 8 x 10 -13 , and clearly outside the scope 
of future experiments. Finally, the e + e~[i + fi~ mode shows a sizeable effect due to the form 
factor that leads to an increase of the width by about 30%. The recent first-time observation 
of this decay mode by the E799 collaboration at Fermilab || furnishes a first experimental 
determination of the branching ratio, namely BR(Kl — > e + e~ / u + /i~)|E xp . = 2.9^2 4 x 10~ 9 . 
Our result is BR(K L — > e + e~// + /i _ )|xheory — (1-30 ± 0.15) x 10~ 9 , consistent with the 
measured value. There is here a strong promise that the reduced error bars resulting from 
the future experimental efforts will permit to show the effect due to the non-trivial form 
factor. This is clearly the most interesting mode for further experimental study. 

As one would have expected, the analysis of the Dalitz decays is enough to pin down 
the predictions for the double lepton pair decays, and the two scenarios resulting from that 
analysis give essentially the same results. 
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